Inelastic electron tunneling spectroscopy at local defects in graphene 
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We address local inelastic scattering from vibrational impurity adsorbed onto graphene and the evolution of 
the local density of electron states near the impurity from weak to strong coupling regime. For weak coupling 
the local electronic structure is distorted by inelastic scattering developing peaks/dips and steps. These features 
should be detectable in the inelastic electron tunneling spectroscopy, d-I/dV^, using local probing techniques. 
Inelastic Friedel oscillations distort the spectral density at energies close to the inelastic mode. In the strong 
coupling limit, a local negative ^-center forms in the atoms surrounding the impurity site. For those atoms, 
the Dirac cone structure is fully destroyed. We further consider the effects of the negative U formation and its 
evolution from weak to strong coupling. The negative U -site effectively acts as local impurity such that sharp 
resonances appear in the local electronic structure. The main resonances are caused by elastic scattering off 
the impurity site, and the features are dressed by the presence of vibrationally activated side resonances. Going 
from weak to strong coupling, changes the local electronic structure from being Dirac cone like including 
midgap states, to a fully destroyed Dirac cone with only the impurity resonances remaining. 

PACS numbers: 73.40.Gk, 73.43.Fj, 03.65.Yz, 68.49.Df 



I. INTRODUCTION 

Graphene has been at the center of attention ever since 
its was first synthesized and studied for its unique physical 
properties While its properties are interesting on its own, 
increasing effort is also being directed towards modifications 
of graphene. Functionalization of graphene has been achieved 
by depositing e.g. H atoms, thus, creating graphane^ which 
is an insulator with a band-gap of the order of 3 — 6 eV. 
Chemical acid treatment may lead to vacancy formation in 
graphene,^ which tends to increase its conductivity due to a 
metallic-like density of electron states (DOS) in the vicinity 
of the vacancies.lSlThe role of single and double vacancies in 
graphene has also been theoretically investigate, showing the 
emergence of midgap states.'^- 

Modifications of electronic states and of the excitation 
spectrum of a given material is crucial for a more effi- 
cient functionalization. Examples of spectroscopies that 
are sensitive to electronic properties are, e.g. photoemis- 
sion and photoabsorption techniques which give access to 
the bulk electronic structure, a nd loc al scanning techniques 
such as atomic force microscopjE^EII and scanning tunneling 
microscop>El(STM). They are employed for studies of spatial 
inhomogeneities^ and local spectral properties."^ 

By studying the response to defects in/on the material im- 
portant spectroscopic information can be accessed.!^ For local 
probes this is a particularly fruitful strategy since it is rela- 
tively easy to move the probe on and off the defect. One thus 
can achieve comparable measurements of the perturbed and 
unperturbed material on one and the same sample. Through 
such an approach effects from potential, charge, andmag- 
netic scattering can be measured from both elasticPSEl and 
inelastic^** point of view. Lately is has become routine to mea- 



sure the inelastic electron tunneling spectrum (lETS) using 
STM. 

In this paper we apply same logic to LETS in graphene. We 
calculate the local density of electron states (LDOS) for elec- 
trons in tight-binding honeycomb lattice which is used as a 
model for graphene. The main results are: 

1 . In the weak coupling limit and using perturbation the- 
ory, the LDOS near the local vibrational impurity ex- 
hibits a kink and logarithmic singularity at the vibra- 
tional mode Oq. The spectral density is significantly 
modified at energies near the vibrational mode. We pre- 
dict those features to be observable in lETS experiment 
using local scanning techniques. 

2. For strong coupling, the atoms surrounding the vibra- 
tional impurity forms negative U centers such that the 
system can be considered as a single impurity problem, 
however, the impurity is effectively spatially extended. 
The LDOS is formed by a series of delta peaks forming 
a single band at negative energy. The result is universal 
in the sense that it is independent of the band structure 
of the conduction electrons, see also She et al. [19] 

3. By coupling the atoms influenced by the vibrational im- 
purity to the surrounding lattice, we study the evolu- 
tion of its LDOS from weak to strong coupling using 
a many-body approach. In the weak coupling regime, 
the Dirac cone is modified by the introduction of elas- 
tic resonances, surrounded by inelastic resonances, sug- 
gesting that the negative f/ -center effectively acts as lo- 
cal impuiity. In the single impurity case, we find that 
all, empty, singly, and doubly occupied, states are pop- 
ulated with a finite probability, which suggests the for- 
mation of local Cooper pair. For increasing coupling. 
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the set of elastic and inelastic peaks move to lower en- 
ergies below the Fermi level leaving a strongly asym- 
metric cone structure around the Fermi level. The Dirac 
cone is eventually fully destroyed in the strong coupling 
limit, leaving two resonances which are broadened by 
the inelastic resonances. 

The paper is organized as follows. First we set up the model 
for the graphene lattice and the vibrational impurity in Sec. |n] 
Then, we move on to discussing the weak coupling limit using 
a T-matrix approach in Sec. [Ill] and the evolution from the 
weak to strong coupling limit using a many-body approach in 
Sec. IV We finally conclude the paper in Sec. |V] 



II. PROBING THE INELASTIC SCATTERING 

We describe the graphene sheet by the nearest neighbor in- 
teraction model 
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where the pseudo-spinor ^,„cr = (fljna ^mcr)' contains the op- 
erators a (b) which annihilate electrons in the A (B) sub-lattice. 

By depositing molecular defect, e.g. CO, on the graphene 
sheet, a local vibrational mode can be introduced. Generi- 
cally, the molecular vibrations cause non-static lattice distor- 
tions. Here, we specifically consider plaquette position of the 
vibrational impurity. An diatomic molecule may, for example, 
be located inside a hexagon in a straight up but slightly tilted 
position.^" The existence of six equivalent positions that the 
molecule can assume, due to the sixfold rotational symmetry 
of the hexagon, may cause molecular rotations, which gener- 
ate local lattice distortions that can be described in term of a 
local bosonic mode coupling to the electronic density at the 
nearest C atoms. Stretching and breathing modes may also 
be envisioned, and especially if the molecule is off-centered 
within the hexagon. Thus, the coupling may be symmetric or 
asymmetric to the near carbon atoms. Here, we shall consider 
both possibilities since the latter can be reduced to effective 
single and double site interactions. 

We, thus introduce (OqB^B, where B^ creates a vibron (local 
bosonic mode) at the energy coq, for the local vibrational mode 
at the position Rq. We describe its coupling to the nearest C 
atoms by 



Xsir,,,) 



where AA/B(r„,) = Aa/b„5(Ro - r„, + Sa/b„), with S„a/b 
are defined in Table |l) whereas Q = B + B' is the vibrational 
displacement operator. Here, the coupling parameters A„ ^ 
Xf,i in general. 

Going over to momentum space via e.g. a„,(j = 
A^^'/'^^ijakcre''* where denotes the number of C atoms 
in the A sublattice, and analogously for the operators on the 
B-sublattice, we can write 



ka 



(2) 



TABLE I: Vectors in momentum space connecting the lattice points. 
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where the potential 0(k) = — f exp (/k • (5„,) such that 
0(k + K±) K, ±Vf feexp{±/(;r/3 — <p)}. Here, the vectors 
are given in Table |l] vp — 3flr/2, tan^ = ky/kx and k = |k|, 
whereas Kj- = ±K = ±2;r(\/3/3, l)/3a. The electron-vibron 
interaction Hamiltonian is in momentum space written as 

^ep = E'J'L^(k,k')'Pk'ae, (3) 

kk'cT 

where A(k,k') = diag{AA(k,k') AB(k,k')} and Xa/b{K^) = 
L« AA/B(r,„) exp [-/(k - k') • r,n]/N. 

III. WEAK COUPLING AND T-MATRIX 

We study the effect of the a weak vibrational impu- 
rity by perturbation theory, which is valid for A^/b,,/? <C 
1. The dressed graphene Green function (GF) G(k,k';z) = 
((^kl^k')) (■^)' suppressing the spin indices, can be calculated 
in terms of the Dyson equation 

G(k,k') =5(k - k')Go(k) + Go(k) ^ S(k, k)G{kX). (4) 
where 

1 I I m\ 



Go(k;z) 



is the bare graphene GF, whereas the self-energy is given by 
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Here, the potential V„,„(z) = /j3"' J^y^l^v - z)G(r„,,r„;zv), 
where we sum over Bosonic frequencies Zv = i2vn/j5, v e Z, 
j3 = l/ksT, and where we have introduced the local Boson GF 
^(2) = {{Q\Q)){z)- In the weak coupling limit, we replace 
both dressed GFs in E by their bare correspondences, using 
Dq{z) — 2(0o/{z^ - COq). Accordingly, the GF is cast in T- 
matrix form in real space 



G(r,r') =Go(r-r') 

+ ^ Go (r - r„, )T (r„, , r„ ) Go (r„ - r') , 



T(r,„,r„) =\d{r,n - r/) - Go(r„, - r,-)V,; 
with the bare real space GF given by 

Go(R 



(7a) 
(7b) 



2n(o I (1) 



(0R\ 



^^(7ocosK-R-i//f'|^ — j 



X I (7v sin 0/; sin K • R + iOy cos cos K • R ) 



(8) 
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Here, H„ '{(o) is the nth Hankel function of the first kind, 
whereas (7,, ; = x,y,z, are PauH matrices and Gq is the iden- 
tity matrix. Here, also R = r — r'. Or — (l)R + n/6, tan^^ = 
(ry — r[.)/{rx — r'^), whereas = Anpv\, with surface den- 



sity p = S/N = k^/An (S graphene area; kc = 2\/2\/l>%l'ia 
large momentum cut off)!22l 

Integration around ±K, yields the retarded potential (with 
obvious notation andx,„„ = p|R„,„|/v'/r) 



VU(«)==7;2A(rm) E / 



m+f{p) 



co — sp — coo + i5 co — sp + coo + i5 
Jo {xmn ) Co COS K • R^n — isJ„ {x,„„ ) [cJ^ sin 9,„„ sin K • R + iOy cos 6,„„ cos K • R] ) pdpX{r„ ) . (9) 



Here, f{x) is the Fermi distribution function whereas «o = 
n{(0o) is the Bose distribution function at coq- 

The electronic structure around the vibrational impurity is 
modified at energies near the inelastic mode ±C0q, where a 
kink and peak/dip is created due to the inelastic scattering off 
the vibrational center Using uniform coupling to the hexagon 
surrounding the vibrational impurity, in Fig. [T] we plot the 
correction to the local density of electron states (LDOS), 
panel (a), and its energy derivative (lETS), corresponding to 
d^I/dV^, panel (b), at Rtip — Rq = a(0,2), for different tem- 
peratures. The LDOS shows non- trivial structure at the vibra- 
tional mode which are more apparent in the lETS as peaks 
around Oo = 15 meV. Similar features are also predicted for 
the case of lETS signatures in d-wave superconductors^^ and 
in topological insulators,'^ as well as for simple metals both 
for vibrationapn and magnetic imputiry.^^ 

The corresponding real space mapping of the lETS is dis- 
played in Fig. [l](c) for energies below, near, and above coq. 
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FIG. 1: (Color online) (a) Change in the local DOS and (b) its energy 
dependent derivative (lETS) at rg = Rq + 2(5i^ (star in the inset) in 
the weak coupling limit, for different temperatures T = 100, 10, 1 
K and vibrational mode (Bo = 15 meV, for uniform coupling to the 
nearest C hexagon, (c) Sequence of lETS maps as function of energy, 
from left to right CO = 12, 15, 17, 20 meV, using T = 10 K and spatial 
broadening F = 2ao/5. We have added an intrinsic broadening of 0.8 
meV in the potential \mn- 



For energies below and above cOq, the presence of the local vi- 
brations generate low contrast, while the contrast grows sub- 
stantially larger for energies around cOq. We expect that the 
presence of the vibrations generates sufficiently large varia- 
tions in the lETS, i.e. d^I/dV^, to be visible in an experimen- 
tal set-up. 

We complete the weak coupling picture by also plotting the 
lETS signatures for asymmetric coupling in Fig. |2j assum- 
ing (a) three, (b) two, and (c) one, C atom being coupled to 
the vibrational impurity. As one may expect, the lETS signal 
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FIG. 2: (Color online) Change in the local lETS at ro in the weak 
coupling limit, for different asymmetric configurations with coupling 
to (a) three, (b) two, and (c) one C atom in the nearest neighbor 
hexagon, and distance from point of measurement, as indicated in 
the upper insets. Left panels show the corresponding lETS maps at 
CO = COq. Parameters as in Fig. [T|(c). 
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is stronger when more C atoms are coupled to the vibrational 
impurity. We also plot different distances between the measur- 
ing point at ro = Rq + 2^i^ and the atom(s) that are coupled 
to the vibrational impurity, clearly showing the oscillatory be- 
havior that is expected due the inelastic Friedel oscillations 
(see insets of the figure, and Fig. [TJ 



IV. EVOLUTION FROM WEAK TO STRONG COUPLING 
REGIME 

We here depart from the T-matrix approximation and con- 
sider the evolution of features from weak to strong coupling, 
i.e. for coupling parameter A^/g/f > 1, using many-body the- 
ory. First, we decouple the Fermionic and Bosonic degrees 
of freedom near the vibrational impurity using the small po- 
laron transformation,^^ that is, constructing the Hamiltonian 
^= e^M'e-^ with 

S=i—Yy]naK^myVmo.P^{-m-B'')- dO) 

We can write the resulting model according to 

+ cooB^B- f^^LA(r„,)^ma) , (11) 



where A(r„,) = A(r,„)/7o^. 

The above expressions are valid for all couplings A^/B(r„,), 
and clearly shows that the presence of the inelastic scattering 
center gives rise to an attractive interaction for the electrons 
residing on the atoms surrounding the vibrational center The 
appearance of the electron-vibron couplings in the first term 
of is due to the fact that S does not commute with a^^bjcr. 



A. Strong coupling limit 

Before we discuss the evolution of the electronic structure 
from the weak to strong coupling regime, we first consider a 
few observations about the strongly coupled system. In the 
strong coupling limit, the system reduces to a single impu- 
rity problem, with the difference to the conventional impu- 
rity problem in that here the impurity is constituted of up to 
six C atoms around the vibrational defect, depending on the 
symmetry /asymmetry of the coupling. For asymmetric cou- 
pling such that the vibrational impurity effectively couples to 
one C atom, the system reduces to a single site problem in 
which the Fermionic ground states can be written, for exam- 
ple, |2) = aj^aj^l), where |) denotes the empty state, assum- 
ing that the vibrational impurity couples to atom « = 1 in the 
A sublattice without loss of generality. The excited states are 



I a) = flifjl) and |0) = |), and the Fermionic energy spectrum 
can be written Ey = — (vA^i)^, v = 0, 1,2. The energy gain 
for doubly occupied site is evident from this result hence we 
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FIG. 3: DOS for the single site problem in the strong coupling/atomic 
limit using cOo = 15 meV and A/Oc = 2 • 10"^ for (a) T = 10 K and 
(b) T = 100 K. Inset in panel (a) shows a close up of the plot in the 
main panel. 



expect this local attraction to play a major role in inducing 
pairing correlations in graphene due to local bosonic mode. 

With the above observations in mind, we write the Hamil- 
tonian of the negative U -site as 

^mp =-Aji^(l +a|^ai<j)a|^aifj. (12) 



In terms of the eigenspectrum of the negative f/ -center, we 
write aia = + oX"^, where X'"' = \p){q\ denotes the 
transition from state \q) to \p) and the factor a = (J^^. We 
can, thus, write 



J^n^p =EoX'^ + £i ^X-^*^ + E2X 



22 



(13) 



The spectrum of the single site is determined through 
the GF G(f,f') = G„^/(f,f')F(f,f'), where G{t,t') ^ 
{~i){T^'„a{t)'i'l^,it')) is the electronic GF and F„(f,f') 



Fnap{t,t') — (=^„a(f)=^„n(f'))vib 



(14) 



is the average over the bosonic degrees of freedom. Here, 

3^na{t) =e"^B^''t^aa{r„)P/cOo^-mB*B,^ a=A,B. (15) 

Following the procedure lined out in e.g. Ref. |26l we calculate 
the generalized function (t = f — t') 

1 

2^ 



(l+2«o) Aa(r„)+AR(r„) 



Fnap{t,t') =exp<! 

-2A„(r„)A^(r„)((l+«o)(l+^^-""') 
+ no(l 

'iving the Fourier transformed GF 

Qr ^^^^ ^^-(l+2«o)[4(r„)+A2(r„)l/2(i)o- 

x£/„(OoK^'**/2G^,,(a)-«t^), (17) 



(16) 
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where (Bq = 2?La{r„)Xp{r„) no{l +no)/<Bo' ^"'^ where I„{x) 
is the modified Bessel function. Thus, for the single site prob- 
lem given by Eq. ([12]), the electronic ground state is in the 
atomic limit given by the GF 



1 — (a\„a 



(0 



5 co + 311, 



(18) 



5 > 0. Setting {a\^aia) = 1, which corresponds to the double 
occupied configuration, we reproduce the analogous spectrum 
found in Ref . J_9. for vibrational impurity on surface of topo- 
logical insulator, i.e. a series of sharp peaks centered around 
the two-Fermion energy — 3A^j . This is shown in Fig. gforCa) 
r = 10 K and (b) T ~ 100 K, also showing that more inelas- 
tic side peaks become activated with increasing temperature, 
as expected. Similar conclusions hold for all our considered 
cases with = 1 , . . . , 6 C atoms coupling to the vibrational 
impurity, with Fermionic ground state consisting of 2N elec- 
trons. 



B. Evolution from weak to strong coupling 

Considering further the single site problem, now in pres- 
ence of the surrounding lattice, we write the transformed lat- 
tice Hamiltonian as 



(19) 



where the coupling between the negative U -center and the lat- 
tice is given by 



^r-L'k(l-e 



(X^Vc7X^'')^ka+^-C., (20) 



^2(7 



with fk = -ri:Li«**'"'^^"V\/^' such that tk±K ~ 
±v/rfee±'(''/^-''')+*'"i /V^V. The negative iZ-center, hence, 
couples to the surrounding lattice with an effective hybridiza- 
tion fk which is renormalized by the momentum P of the local 
bosonic mode. 

We capture the evolution from the weak to strong cou- 
pling limit by solving the equation of motion for the many- 
body operator GF = {{X"\X^)){z), for the transitions 
a,b = (Off), ((t2), self-consistently in mean-field approxima- 
tion under the self-consistency condition that the occupation 
numbers A^o + Ln- +N2 — 1 . The occupation numbers are 
calculated using^^ 

^o--^ImE/[l-/(®)]G(Wo(»M», 

Nc,^- ^Im I [/(a))G^^^o(®) + [1 - f{(oW^22oi(o)]dco, 

(21b) 

N2^--lmY^ /"g^ (21c) 

Due to the inherent spin-degeneracy and absence of a coupling 
between the spin-channels, the GF reduces to a 2 x 2-matrix 



equation. To second order in fk and X, the result is given in 
terms of the retarded GF 



(to) = CO- A-PE(a))(l + (jO 



(22) 



where A = diag{Ai A2}, A„ = E„-E„^u P = diag{Pi P2}, 
Pi = No +N1/2, P2 = N1/2+N2, Ni = l^aNa, whereas the 
self-energy r(co) = Ln=i,2^''"H'^) is given by 



S(i)(co) = 

S(2)((0) = 



2to 



2 log 



insignco 



(23a) 



-47t 



/(co) CO CO + X^ 



to sign CO. (23b) 



The contribution E^^' account for fluctuations on and off the 
negative U -center, essentially caused by the presence of the 
surrounding lattice, showing a cubic correction to the LDOS. 
The second contribution, is generated by fluctuations on 
and off the negative [/-center due to the coupling between the 
local vibrational mode and the Fermionic degrees of freedom. 



Equations (23 1 should be solved self-consistently, however, 
we can make a few observations on the expected behavior of 
the electronic structure. For weak coupling, the bare excita- 
tions Ey = — (vAai)^ 0. Thus, for low energies, such that 
co/Dc ^ 1, we can neglect the self-energy E'^^ and approxi- 
mate the first self-energy by E^'^ « —2a). Then, the denomi- 
nator of Gi^ is given by 

2 2 

3to2-2co^ A„ + ]~[A„ =3(co-A+)(ffl-A_), (24) 

n=\ n=\ 

where A± = —(4^ /3. This is found by observing 

that Ai = — and A2 = —31^, such that L„ A„ —Aa^ and 
Hn^n = 3A^. Here, we have, moreover, used that P\+P2 = 
No + Y,aN(y + N2 = 1, by charge conservation, along with 
Pn- 1/2, c.f Fig.g 

As the coupling is increased, the non-linear components in 
the self-energies play an increasingly important role for the 
positions of the poles, such that we cannot any longer make 
use of Eq.d24]). 

In Fig. |4|(a) — (d) we plot the evolution of the LDOS on 
the negative fZ-center from weak to strong coupling regime 
for low (bold/black) and high (faint/red) temperatures. In the 
weakly coupled system, panels (a), there are two main (elas- 
tic) peaks near the Fermi level, corresponding to A±, c.f. Eq. 
(21a) (24 1. For low temperatures there is tiny signature of a vibra- 

15 meV. For higher 



tional side peak at about (O — —(Oq 
temperatures, these vibrational signatures become more ap- 
parent, as one should expect since those modes are thermally 
activated. 

For increasing coupling the main elastic features remain, 
however, shifted to lower energies. They become increasingly 
broadened since the level width is cubic function of the en- 



ergy, c.f. Eq. (23a I. Moreover, the presence of the vibrational 
side peaks also become more visible in the LDOS, even for 
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FIG. 4: (Color online) Evolution of the LDOS at the negative U- 
center from weak to strong coupling regime. Here, X/Dc = {5 • 
10""*, 1 • 10"^ 5 • 10"^ 1 • 10^2}, CB() = 15 meV, and T = 10 K 
(bold/black) and T = 100 K (faint/red). The inset in panel (c) shows 
the full LDOS at T = 100 K. 



low temperatures. In both cases illustrated by panels (a) and 
(b), the coupling is weak enough to preserve the overall Dirac 
cone, apart from the presence of the resonances. 

For even stronger coupling, panels (c) and (d), the Dirac 
cone is fully destroyed and only the peak features, caused by 
the elastic and inelastic scattering, remain. Finally, in the 
strong coupling limit, panel (d), there only appears a dou- 
ble peak structure, where the peaks coiTespond to the singly 
and doubly occupied states. For high temperatures, the vibra- 
tional side peaks effectively act as a thermal broadening of the 
main peaks. The discrepancy with the situation illustrated in 
Fig. [3]can be understood from the fact that we here take into 
account fluctuations to both the singly and doubly occupied 
states, hence, there is a finite likelihood that even the singly 
occupied state becomes populated. This is typical feature of 
any many-body description, and it emphasizes the fact that the 
charge is partially distiibuted among the available states. 

We finally comment on the evolution of the Fermionic state 
of the negative fZ-center from weak to strong coupling regime, 
represented in terms of the populations numbers A^„, n = 
0,1,2, c.f. Fig. |5] In the weakly coupled system, the energy of 
the single electron fluctuations A„ — —{nXf- + (n — 1)^/1^ = 
— (2n — 1 ) Ues below but close to the Fermi level, c.f. Fig. [4] 
(a), (b), such that the system is open for fluctuations between 
the (four) states. This property is verified by the occupation 
numbers, in that all A^„, n = 0, 1 , 2 are finite. This suggests oc- 
currence of local Cooper pair formation near the vibrational 
impurity, which will be the topic of a future publication. 

In the strongly coupled limit, on the other hand, the set of 
elastic and inelastic transition energies are far below the Fermi 
level, c.f. Fig. |4](d), such that the the population number 
A^o approaches zero. In other words, the negative t/-center 
acquires a Fermionic ground state which is a mixture of the 
singly and doubly occupied states. The coupling between the 
negative [/-center and the surrounding lattice, thus, generates 
a more intricate electronic structure than what is suggested 
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FIG. 5: (Color online) Evolution of the occupation numbers A'o (tri- 
angles), N\ (pentagrams), and Ni (bullets), at the negative IJ -center 
from weak to strong coupling regime. Here, coq = 15 meV and 
r = 10 K. 



by the atomic limit physics where the negative [/-center is 
decoupled from the lattice. 

In the intermediate regime, there is a cross-over regime, or 
possibly a phase transition, c.f. crossing of population num- 
bers near XjDc > 10^^ in Fig. [s] where the occupation num- 
bers of the empty and doubly occupied states evolve monoton- 
ically decreasing and increasing, respectively, with the cou- 
pling strength A, whereas the single Fermi on state(s) remain 
constant. 



V. CONCLUSIONS 

We have theoretically studied the effects of vibrational im- 
purity adsorbed onto graphene, specifically the inelastic scat- 
tering properties. We find in the weak coupling regime, that 
the perturbed LDOS in the vicinity of the vibrational impurity 
acquires peaks/dips and steps at the energy of the vibrational 
mode. The spectral density distortions around the vibrational 
mode is spatially extended showing inelastic Friedel oscilla- 
tions, in analogy with the findings for surfaces of metallic 
materiald^iJ21IMl and topological insulator.'^ 

By employing a many-body approach, we study the evo- 
lution from weak to strong coupling regime. In the weak 
coupling regime, an elastic mid-gap resonance emerge, sur- 
rounded by inelastic side resonances, at half the energy of 
the single electron fluctuations between the negative U -center 
and the surrounding lattice. Finite occupation of all Fermionic 
states, empty, singly, and doubly occupied state, on the neg- 
ative [/-site, near the vibrational impurity in the weakly cou- 
pled system, suggests local Cooper pair formation. The as- 
pects of this physics will be the topic of a future publication. 

For intermediate coupling strength the peak structure is 
severely distorted and pushed below the Fermi level, leav- 
ing a strongly asymmetric Dirac cone around the Fermi level. 
The Dirac cone is eventually destroyed in the strongly cou- 
pled regime, in which the electronic structure acquires a band 
formed by the collection of elastic and inelastic resonances. 

We believe that our findings should be within the scope of 
present experimental local probing abilities using e.g. STM 
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or atomic force microscopy. 
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